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Abstract 

Relativistic massive bosons with spin one are considered in several quantization 
schemes. In all of them the system is shown described by a non-Hermitian Hamil- 
tonian H ^ and helicity operator A. Constructively we show that in all of 
the contemplated schemes both these operators H and A prove simultaneously 
VT— symmetric, i.e., pseudo-Hermit ian with respect to a certain not too compli- 
cated indefinite pseudo-metric operator V = V'. 
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1 Introduction 



The majority of textbooks on Quantum Mechanics illustrates the performance of 
the formalism via a non-relativistic point particle of mass m which moves in a one- 
dimensional potential V(x). On the basis of the principle of correspondence the time 
evolution of such a system is quite often easily described as generated by the ordinary 
differential self-adjoint Hamiltonian operator H = p 2 + V(x) where p 2 = —d 2 /dx 2 
in units h = 2m = 1. Much less attention is usually being paid to the possibilities 
of the transition to relativistic kinematics [1], i.e., to the Klein-Gordon equation for 
the spinless bosons, to the Dirac equation for the simplest fermions, to the Maxwell 
equations for the massless photons, to the Proca equation for the massive bosons 
with spin one [2] etc. 

It is well known that the correct description of relativistic systems requires the 
use of the full-fledged formalism of relativistic Quantum Field Theory. Another, 
"hidden" reason of the overall reluctance of working with relativistic kinematics 
even in an approximative regime of Quantum Mechanics is that it often requires a 
transition to non-self-adjoint generators of the time evolution, i.e., to the Feshbach- 

Villars Hamiltonian operator i/( FV ) ^ (ii"( FV )) in the spinless case [3] etc. Appar- 
ently, the naive use of the similar Hamiltonians may lead to inconsistencies. For 
example, negative probabilities of the localization are traditionally being mentioned 
as occurring for the relativistic spinless boson at a point x [4], etc. 

The formalism based on the FV-like Hamiltonians f/~( FV ) need not necessar- 
ily be mathematically inconsistent [5, 6]. In fact, the use of H = H^ FY "> may be 
re-interpreted as one of the most characteristic applications of the so called PT- 
symmetric version of Quantum Mechanics as originally proposed by Bender and 
Boettcher [7] . Thanks to a concentrated effort and subsequent debate (cf . its sample 
in [8]-[14]) it became clear that even the manifestly non-Hermitian PT-symmetric 
operators of observables O (characterized, for our present purposes, by their prop- 
erty 

= vov~ 1 ^o (1) 

with a suitable "indefinite metric" operator V or with some of its non-metric alter- 
natives [15]) may fit in the overall scheme of Quantum Mechanics. 

On this formal background and in the language set already by Scholtz et al 
in early nineties [16] it has been recognized that in the relativistic spinless case the 
Hamiltonian if ( FV ) may be treated as self-adjoint with respect to some less standard 
inner product in Hilbert space. In the other words, the use of the non-Hermitian 
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observables may often be made fully acceptable via an introduction of a certain 
nontrivial, positive definite "physical" metric operator ^ I [12] - [16]. 

All the "physical" choices of © must make all the underlying observables (say, 
O = ij( FV ) etc) self-adjoint in some new, ©—dependent Hilbert space, 

0t = q @-i ^ 0> = t > 0. (2) 

In this sense, no modification of the postulates of Quantum Mechanics is needed. 
Nevertheless, in order to avoid confusion the authors of the review [16] recommended 
that one should speak about quasi- Hermiticity of O whenever ^ I in (2). 

For the relativistic spinless system the "weakened" Hermiticity property (2) of- 
fered a natural explanation of the negative-probability Klein-Gordon puzzle [6]. In 
what follows we intend to discuss some aspects of the transition to the "next" case 
of the massive field with spin one. The main emphasis of our paper will be put on 
the existence of the VT— symmetry features of this system and on the emergence of 
their specific ambiguities which arise in connection with the presence of a nontrivial 
spin s > which mediates a coupling of channels [13]. 

Our present results will be separated into their classical part (section 2), quan- 
tization part (section 3) and VT — symmetry-related part (section 4). Our observa- 
tions starting from the classical case may be perceived as separating the Lagrangian 
and Hamiltonian formulations, with the core of our attention concentrated upon 
the latter language. Under both the classical and relativistic kinematics our study 
will be sub-separated into the economical "non-Dirac approach" without any con- 
straints (cf. subsections 2.3 and 3.1 or 4.1, respectively) and the more common 
"Dirac approach" using constraints (cf. the respective subsections 2.4 and 3.2 or 
4.2). 

In the latter, Dirac-theory context there still exists certain freedom in the choice 
of the method of the quantization. Three possibilities, viz., "version A" (cf. 2.4.1 
or 3.2), "version B" (cf. 2.4.2 or 3.3) and "version C" (cf. 2.4.3 or 3.4) are studied 
in more detail in what follows. 

As long as the key purpose of our study lies in an extension of the concept of 
VT— symmetry (1) to the equations controlling the motion of the relativistic vector 
bosons, section 4 represents in fact a climax of our considerations. A number of 
explicit formulae is derived and displayed there including the matrix forms of the 
Hamiltonian as well as of the related pseudo-metric operator V. 

Our last section 5 brings a brief summary of our results while Appendix A adds 
a brief discussion of a few parallels between our present (massive) model and its 
more common massless (i.e., electromagnetic) counterpart. 
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2 Bosons with spin one 



Our project has been inspired by the description of the massive vector boson system 
by Taketani, Sakata and Tamm [17] as summarized by Nikitin et al [18]. Our sample 
of relevant references shouldn't also omit Kemmer [19] (in the context of the so-called 
Kemmer-Duffin-Petiau equation) or Shay and Good and Weinberg [20] (in a slightly 
different context) and many others, with a nice summary of the field offered by 
Labente [21]. 

2.1 Lagrangians 

Let us start now our present discussion from the classical case and assume that the 
dynamics of the classical complex vector field A of mass m is given by its Lagrangian 
density 

C = m 2 A^ - \ (d^A v - cU M ) {dPA v - d v A») . 

The related action S — J Cd A x has to be stationary 5S = 0. With the current 
abbreviation □ = —g flu d Mi d u where g^ v = diag(l, — 1, — 1, — 1) (and in units such 
that c = 1 and x = t) this gives the usual Euler-Lagrange equation 

d il F tlv + m 2 A" = 0, = &*A V - d v A» (3) 

plus its complex conjugate. Its differentiation and summation reveals that the indi- 
vidual components of the field have to fulfill Klein-Gordon equation subject to the 
free Klein- Gordon- type constraint, 

(m 2 -ri)A li = 0, d v A v = 0. (4) 

Consequently, only three of the four components of the field are linearly independent. 
The equation itself is obviously satisfied by the wave solutions distinguished by the 
sign a — ± of the energy, 

A^x) = N a (p)u^p)e-^, Po = au, cu = |^P 2 + m 2 |. (5) 

An additional constraint acquires the form pu a (p) = 0. For the sake of definiteness 
the three linearly independent four- vectors u a {p) may be chosen in the following 
form 

u aA = (0, u ajl ) T , u a)2 = (0, u aj2 ) T , u Q;3 = (a^j-, j^jl • (6) 
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The space-like components (i.e., polarization vectors) in the first two items are 
perpendicular to p while the third one is chosen parallel to p. We may fix 



, x (o, -p 3 ,P2) T ( s (p| +ph-pm,-pm) T 

U ",HP) = / = U a,2(P) " 



pi 



Alternatively, the operator of helicity 



h(p) = ^ = nS 

IpI 



\J{pI +pi) 2 +P1P2 +pIpI 



I — n 3 n 2 \ 

n 3 —rii 
\ -n 2 Tii / 



(7) 



(8) 



becomes defined in terms of the three generators of the three-dimensional represen- 
tation of the rotation group, 
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Its action on the vectors (6), [hu a ^ = iu aj2 , hu a 2 = —iu a ^ and hu a;3 = (with 
subscript a — ±)] prompts an introduction of the helicity-numbered eigenvectors 
u«(p,±l), 

hu a (p, ±1) = ±u Q (p, ±1), hu a (p, 0) = 0, (9) 

with an explicit representation u±(p, +1) = (w±,i + iu± t2 ) / \/2, u±(p, —1) = (ii±,i — 
iu±,2)/V2 and u±(p, 0) = w±,3(p). Combining (5), (6) and (7) we get 



i.e., waves traveling at a fixed energy and helicity. 



(10) 



2.2 Hamiltonians 

There exist, basically, two ways of the derivation of the description of the time 
evolution of the vector boson fields from the first principles. In one of them one 
eliminates all the (redundant) degrees of freedom from the very beginning. We shall 
call this approach "economical" because the number of its field "parameters" is, in 
some sense, minimal and equal to the number of the degrees of freedom of the field 
(cf. also subsection 2.3 for more details). 
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An alternative recipe (cf. the subsequent subsection 2.4) carries the name of 
Dirac [22]. It employs the language of constraints and its key merit may be seen in 
its more formal, "algorithmic" character. It offers a virtually unambiguous recipe for 
a systematic treatment of the field in question. Both its physical and non-physical 
components are treated "on an equal footing" . 

The key purpose of our present paper lies in the discussion of the latter two 
approaches. We intend to illustrate (by construction!) that both of them remain 
equivalent on the classical level while a number of open question can arise during 
their quantization. 

Before we fully concentrate upon the questions of specific differences between 
quantization recipes for vector bosons (in which the classical fields should be replaced 
by operators) we intend to emphasize that several important technicalities must 
already be re-analyzed in a preparatory step, on the classical level. On this level the 
first relevant aspect of the problem may be determined immediately as lying in the 
existence of the "economical" possibility which looks nonstandard and which seems 
to be in an apparent conflict with the more modern formulations of the problem. 

In the modern (usually called Dirac's) theory the quantization is treated with 
the BRST and Batalin-Vilkovisky methods [23]. For a more detailed elucidation 
of the topic the reader might consult the related literature where the quantization 
of the relativistic (free) particles is reminiscent of Feynman's fifth time formalism 
and where the theory involves the quantization of the reparametrization-invariant 
models on the world line, etc. 

The latter approach leads to the Klein-Gordon equation emerging as a constraint 
(analogous to the Gauss law in QED) and it opens new perspectives, i.a., in super- 
symmetry. At the same time it does not represent the only possibility and it may 
often be complemented by the less subtle approach where the Klein-Gordon equa- 
tion results from the "economical" process where one starts from an elimination of 
all the redundant degrees of freedom. Thus, once we succeed in a minimization of 
the number of the fields we are not forced to impose any auxiliary constraints. We 
shall add more details in paragraph 2.3. 

One could stress that the minimization strategy circumvents (or at least clarifies) 
certain shortcomings of the use of the Dirac's eight-dimensional Hamiltonians by 
their replacement by their compact, six-dimensional "economical" alternatives. Vice 
versa, the specific merits of the eight-dimensional formulae will be recollected in 
paragraph 2.4. We shall clarify there the distinction between several alternative ways 
of the detailed application of the Dirac's recipe. We shall distinguish between certain 
"approach A" (cf. paragraph 2.4.1 for more details), "approach B" (cf. paragraph 
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2.4.2) and "approach C" (cf. paragraph 2.4.3). We believe that in this way an 
enhanced transparency will be assigned not only to the essence of the (second-class) 
constraints but also to the subsequent questions of the quantization. 

2.3 "Economical" formulation 

Using the standard recipe and taking and A^ as generalized coordinates we get 
the generalized momenta using the standard recipe, 

dC j j dC 

d(d A m ) d{d A ) 

dC - dC 

Um= a \ = d o A m ~ d m A , n = - = 0. (11) 

We see that the time-coordinate momenta Ilo and llo can be omitted whereas their 
respective field conjugates can be eliminated using the equations of motion (3), 

A = — 2 (dJL), A = \ (d rn U m ). (12) 

After a straightforward calculation the density of the Hamiltonian of our classical 
vector field can be found, 

H = U m (d A rn ) + n m (d A rn ) -£ = 

= n m n m + Ar(<9 m n m )(<9 n n n ) + m 2 A m A m + (d m A n )(d m A n ) - (d m A m )(d n A n ). 



m 



Classically, Euler-Lagrange equations (4) correspond to the dynamical equations 

A m = {H,A m } P , n m = {ij,n m }p (13) 

where {., .}p denotes the Poisson bracket and where H — J Hdx. 



2.4 Dirac theory using constraints 

We already mentioned that there are several possibilities of the strictly formal tran- 
sition from Lagrangians to Hamiltonians. Sometimes, one has to interpret the im- 
possibility of the elimination of the generalized velocities as an emergence of certain 
phase-space constraints (to be denoted as 6). They are called "primary constraints" 
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and will be marked by the superscript ® . We shall have to deal with a pair of them 
in our vector-boson model. 

In addition, the natural requirement of the time-independence of the constraints 
9^ leads, in principle, to the emergence of the further series of constraints denoted 
as They are called "secondary" and we will have a pair of them in what follows. 

In version "A" of the formalism (cf. paragraph 2.4.1 below), the Hamiltonians 
will be constructed, using the technique of the Lagrange multipliers, in a way which 
employs the primary constraints 9^ only. In version "B" of paragraph 2.4.2, the 
knowledge of all the constraints 9^ and 9^ will be employed giving an alternative 
family of the Hamiltonian functions. Finally, the approach "C" of paragraph 2.4.3 
will re-define the brackets and enable us to prepare and perform the quantization in 
an internally most consistent and compact manner. 

2.4.1 Approach A 

In the approach "A" one has to apply the Lagrange multipliers in order to respect 
the existence of the two primary constraints. These constraints read 

0p) = n = o, = n = o (14) 

and emerge in the process of the transformation of the Lagrangian into the Hamil- 
tonian (cf. also the formulae of eq. (11)). An "extended" Hamiltonian will be 
then sought as containing a pair of the Lagrange multipliers A® and A®, 

+ \®9® + \®W. (15) 

The symbol TC° stands for the standard, non-extended Hamiltonian function, 

H° = U m (d A m )+Tl m (d A m )-C= (16) 
= II m n m + m 2 A m A m - m 2 A A + (d m A n ) (d m A n ) - 
-{d m A m ){d n A n ) + U m {d m A ) + fl m {d m A ). 

In the next step one demands the time independence of the primary constraints, 

{H^ A \ 0«} P = 0, {H^ A \ 9^} P = 0, (17) 

where 

= J dx7i {K) . (18) 
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In our present model such a requirement gives rise to the mere doublet of secondary 
constraints, 



0(°) = -m 2 A - (d m U m ) = 0, = -m 2 A - (d m U m ). 



(19) 



In the subsequent step we have to demand the time independence of the secondary 
constraints, too. In a way paralleling eq. (17) we obtain the Lagrange multipliers 
A® and A® as a solution. In this way we arrive at the final Hamiltonian 



The process of construction guarantees that all the equations of motion respect all 
the primary and secondary constraints. The same information is of course carried 
also by equations (12) which were employed in the "economical" approach of para- 
graph 2.3. 

2.4.2 Version B 

In the previous approach A, the secondary constraints were clearly discriminated 
in comparison with the primary ones. Fortunately, in a way outlined, e.g., in the 
popular textbook [24] we can also introduce all the constraints in the Hamiltonian 
simultaneously, regardless of their origin. 

First steps of such a slightly more sophisticated approach to the construction 
of Hamiltonians are exactly the same as in the approach A. The differences appear 
after one obtains all the secondary constraints. Having collected all the constraints 
6® and 9^ one decides to employ them all in the recipe using Lagrange multipliers. 
This means that our new postulated Hamiltonian Ti^ will read 



where the symbols A* 11 ) and A^ 1 ) denote two new Lagrange multipliers attached to 
the secondary constraints. All the values of the lambda-multipliers emerging in 
eq. (21) may be again extracted from the requirement of the time independence, 




(20) 



^(B) = n° + \®e® + \®e® + A< n >0< n > + A (II ¥ n) 



(21) 






The resulting Hamiltonian 



H (B) =n o + U ^ dmAm) + fi (d m A m )+ 
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+— 

m z 



m 2 A + (<9 m II m ) m 2 A + (<9 m II m ) 



(23) 



will prescribe the motion equivalent to the one obtained within approach A, on the 
classical level at least. 



2.4.3 Approach C 



A key shortcoming of both the above approaches A and B (which lead to the same 
equations of motion as obtained in the "economical" setting) is that they prove 
less convenient for quantization purposes. Apparently, this would make the Dirac's 
theory manifestly disqualified in comparison with the economical approach. Fortu- 
nately, it is possible (as well as not too difficult) to re-adapt the Dirac's procedure 
for quantization purposes. 

In the first step we calculate the Poisson brackets for all the constraints and 
define a certain auxiliary matrix Ai, 



(M (x, y 
where we renumbered 



J,K 



{0^(x),9^ K \y)} P , J, K= 1,2, 3, 4 



0« = e®, 

0(3)= 0(1) , 



0(2) = 0(H) ? 
0(4) = 0(11) _ 



We find that 



M JK (x,y) = M JK 5^\x-y) 
with the boldface symbol M denoting an antidiagonal matrix, 

/ m 2 \ 



(24) 



(25) 



(26) 



M 



V 



o 




2 

-TO 






2 

TO 





TO 

-m 2 



/ 



M 



m 4 



(27) 



As long as the matrix M is non-singular and invertible we can easily insert 

Mj 1 K ^,y) = Mj 1 K 6^(x-y) (28) 
in an explicit definition of the Lagrange multipliers of approach B and obtain 



(29) 



K=l 
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The numbering of As is inherited from 8s (cf. eq. (25)) and gives the final compact 
formula for the Hamiltonian, 

H (B) =n °- { H °, fl (Jf) }pM^ (J) , H°= [ dxH°. (30) 

J,K=1 J 

Still, our work is not finished yet since the consistency of the recipe still requires a 
replacement of the Poisson brackets by the so called Dirac brackets, 

{X,Y} D = {X,Y} P - £ ( dwd^X,d {J \w)} F Mf K {w,z){d {K \z),Y} F . (31) 

J,K=1 J 

Fortunately, one can show very easily that the equations of motion which result from 
approach B are practically equivalent to those which can be obtained from eq. (31). 
Roughly speaking, there emerge differences reducible to certain combination of terms 
which would vanish due to the constraints. Ref. [24] may be consulted for details. 

The main advantage of the replacement of the Poisson brackets by their Dirac 
alternative is that the Dirac brackets annihilate all the constraints, 

{X,6 {J) } D = for Je {1,2,3,4} and VX. (32) 

Such a property will prove extremely useful during quantization since the use of the 
Dirac brackets (in place of the Poisson ones) enables us to work with the "core" 
Hamiltonian 7i° only. 

3 Quantization in Heisenberg picture 
3.1 "Economical" formulation 

In the Heisenberg picture the quantization of the system yields a replacement of 
the generalized coordinates A m and Yi n and of their complex conjugates by the 
non-commutative operators which satisfy the canonical equal-time commutation re- 
lations, 

[4n(x), n n (y)] = iS nm 5 (3 \x - y), 
[An(x), H n (y)] = i6 nm 5®(x-y), n,m = 1,2, 3. (33) 
They are — up to the factor % — parallel to the Poisson-bracket rules. 
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The time-evolution of the operators is determined by the modified equations 
(13), 

A k = -i[H,A k ] = -U k + — d k (d n U n ), 

m z 

U k = -i[H, n fe ] = m 2 A k - AA k + d k {d n A n ), 

A k = -i[H,A k ] = -U k + —d k (d n U n ), 

m z 

rife = -i[H, U k ] = m 2 A k - AA k + d k {d n A n ). (34) 

An arrangement of A and n in a single six-component set \1/ T = (mA, ill) is some- 
times being used [25]. 

3.2 Dirac theory — approach A 

In place of the six-dimensional Hamiltonians H pertaining to the preceding approach 
let us now switch to the so called Dirac's approach [22, 24] where, as we already 
outlined, certain more-dimensional Hamiltonians with redundant components are 
being used. 

The form of the commutators should result from their Poisson-bracket prede- 
cessors. Nevertheless, we have to keep in mind the existence of the additional, 
"redundant" degrees of freedom now. Thus, an additional pair of operators with the 
non-trivial commutators emerges, 

[A (x), n (y)]=^(x-y), 

[A(x), n (y)]=^ (3) (x-y). (35) 
Performing all the calculations we may collect the extended set of equations, 

A k = -i[H^\A k \ = -fife - (d m Ao), 

rife = -i[H^ K \ lb,] = m 2 A k - AA k + d k (d n A n ) - (<9 m rr ), 

A = - i [H^ A \A ] = -(d m A m ), 

n = -t[H^\ n ] = -m 2 A - (d m Tl m ). (36) 

For the sake of brevity we do not display the complementary, complex conjugate 
version of these equations. 
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3.3 Dirac theory — version B 

The Hamiltonian enters formulae 

A k = -i[HP\A k ] = -U k + — 2 d k {d n tl n ) + (d m A ), 

ri* = -i[H (B \ U k ] = m 2 A k - AA k + d k {d n A n ) - (9 m n ), 

A = -t[H^\A ] = -(d m A m ), 

rio = ^ n ] = +m 2 A + (9 m n m ). (37) 

In essence, the practical consequences of the procedure remain virtually the same 
as before. 

3.4 Dirac theory — approach C 

The major change comes with the necessity of using the "Dirac-originating" com- 
mutators in place of the "Poisson-originating" ones. This leads to the following two 
nontrivial rules 

L4 (x), n (y)] = 0, 

L4 (x), B m (y)} = -^5 (3) (x - y) (38) 

(cf. relations of eqs. (33) and (35)). In the light of the latter "distortion" of the 
commutators we are allowed to employ the Hamiltonian H° and to arrive at the set 
of the equations 

A k = —i[H^ B \ A k ] = -fi fc - \d k {d n U n ) - 2(d m A ), 

ri fe = -%\E^\ life] = m 2 A k - AA k + d k (d n A n ) 

A = -t[H^ B \A ] = -(d m A m ), 

tl = -t[H^ B \U } = 0. (39) 

We may observe and conclude that all the available recipes are equivalent because 
they differ only by the terms which are proportional to certain combinations of the 
constraints. 
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4 Time-evolution and VT— symmetry 



4.1 Six-dimensional Hamiltonian 

After the "economical" quantization we arrive at the formal Schrodinger equation 

<^|*> = J ff|*> (40) 

which prescribes the time-evolution of the operators in our system. The quantum 
Hamiltonian H = ij( VB ) of the vector boson is a matrix 



-m + 



grad div 



+ o "' I (41) 
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-d 3 d 1 


-<9 3 <9 2 


-9 2 - c 2 











/ 




(cf. [25]). In the coordinate representation it has the 6 x 6— dimensional and mani- 
festly non-Hermitian partitioned matrix structure, 



H 



As a V— pseudo-Hermit ian operator, 

fl-t = -pHV~ l ^ H, V = 

it allows a quantization in a finite box of a side-length L. Using the periodic bound- 
ary conditions imposed at its walls one arrives at the discrete bound-state energies 
on the mass shell, 

Po,n = iy/p 2 +m 2 , p jta = — n 

with n being a vector of integer components. Vice versa, the continuity of the 
energies can be restored in the L — > oo limit whenever necessary. 

At a finite L the six independent eigenvectors with eigenvalues ±p may be 
most easily constructed in terms of the above-introduced functions in momentum 
representation, 

i) - I mA^ipJ), k = 1,2,3 _ 
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Here we abbreviated 

, 4<± , (p , 1) = 4%H^> M 

A«(p,3)=4 ±) W 

with three linearly independent auxiliary amplitudes A ( f 1 \p) introduced originally 
as perpendicular to the momentum four-vector, 

4 ±} ( P ) = (o, ei ), 4 ±} (p) = (0,62), 4 ±} = (— ■ ( 42 ) 

\Po |P|/ 

The unit vectors e 1; e 2 are chosen as mutually perpendicular, (ei,e 2 ) = 0, and 
perpendicular also to the momentum vector p, 

A (±) ( P ) = 4 ±} (p) ei + 4 ±} (p)e 2 + 4 ±} (p)p • 

In order to make our solutions (confined in a box!) less ambiguous we introduce 
another observable which commutes with H, 

A i,j\P) - A i+3,j+3{P) - —Z H , 

z p 

A w (p)=A l+3j (p) = 0, i,j = 1,2,3. (43) 

This is a block diagonal operator which distinguishes the spin projection into mo- 
mentum direction. In terms of physics it corresponds to the helicity of our vector 
bosons. It is easy to verify that the latter operator is also pseudo-Hermitian with 
respect to V, 

At = VAV- 1 . 

We obtain a complete right-action Schrodinger equation for the fixed-helicity states, 

H\ ± cu, h) — ±u;| ± cu, h), A\cu, h) — h\cu, h), h = ±l, 
as well as the similar relation for the action of our observables to the left, 

({±u,h\H = ±u;{{±u},h\, ((±tu,h\A = h((±u,h\, h = ±1, 0. 
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In a way described more thoroughly in ref. [13] the pseudo-Hermicity of both our 
observables can immediately be used to express the double-ket-marked eigenvec- 
tors of the conjugate-operator pair of W and A* in terms of the single- ket-marked 
eigenvectors of the pair of operators H and A, 

\±u,h)) = P\±u,h)Q P0 , h . (44) 

One can normalize the vectors \po,h) and \p' ,h')) in such a way that they form a 
biorthonormal set, 

<(Po, h\p' Q , ti) = 5(p Q - p' )5 hh > 
(cf. ref. [13] for more details). 



4.2 Alternative eight-dimensional Hamiltonians 



In the Dirac case we can write Hamiltonians in the partitioned matrix form, too. 
All such matrices will be eight-dimensional. 

In approach A, the matrix form of the Hamiltonian is compatible with the related 
equations of motion (36), 



( 

_ m _l_ Yl — § rad div 

m m 



-m + 



V 



div 




2 grad div 
m 




div 



grad \ 

grad 



—m ) 



(45) 



It is amusing to notice that the approach B will give the matrix Hamiltonian in a 
slightly different form, 



H (A) 



( 

— m _|_ Yl _ e rad div 



V 



div 




— m 



div 



grad \ 

grad 



m J 



(46) 



(cf. equations of motion (37)). Still another version of the Hamiltonian will emerge 
from the approach C, 







-m + ^ 

m 



grad div 

TO 



div 




-m - -2 grad \ 





) 



(47) 
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(cf. equations of motion (39)). All the latter matrices are non-Hermitian. However, 
they are also VT— symmetric in the way which employs the same, constant and 
partitioned matrix operator 



V = 



( o h o \ 

h o 

-1 

V o o-i o / 



(48) 



in all the three contexts. It is worth emphasizing that the use of the Hamiltonian 
pertaining to the Dirac approach C seems definitely preferable from the point of 
view of quantization. 



Once we assume that the pair of a quantum Hamiltonian H and of some auxiliary 
"observable" A are independent and diagonalizable non-Hermitian operators with 
discrete spectra, we may study this pair as a special mathematical realization of 
a VT— symmetric quantum model [13]. In our present paper we paid particular 
attention to the related necessity of the explicit knowledge of a sufficiently simple 
auxiliary generalized parity operator V . In this sense we have shown here that the 
classical Proca's field offers a new and feasible example of a VT— symmetric system, 
this time with an immediate physical vector-boson interpretation. 

5.1 Proca system with pseudo-Hermit ian observables 

Let us summarize our preceding text as a constructive presentation of the triplet of 
operators H, A and V on a Hilbert space V which fulfills the relations 



where V is an extremely simple though nontrivial and positively indefinite Hermitian 
automorphism on V. 

In the related literature, the PT— symmetric operators H and A are also often 
called pseudo-Hermitian or V— pseudo-Hermitian. One of the most important con- 
sequences of their VT — symmetry is well known to lie in the facilitated possibility of 
the construction of the physical or "representation" Hilbert space and metric using 
the factorization ansatz O = CP where C may be called a "charge" [12]. 



5 Summary 



ift = VHV' 1 , A f = VAV~ 



i 



(49) 
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Let us repeat that having assumed the simultaneous diagonalizability and & 
non- degeneracy of the (real or complex) spectra of both H and A we can, at least in 
principle, find their mutually bi-orthogonal left and right eigenvectors, 



H\n, h) = E n \n, h), 



((m,h\H = {{m,h\E, 



(50) 



A|n, h) — h\n, h), ((m, h\A = ((m, h\h, 
as well as their formal spectral representation, 



(51) 



H 



{{n, h\n, h) 



((n,h\, 



(52) 



n 



Due to (49), the left and right eigenvectors are mutually related by V 




(53) 



V\n,h) ~ p nth \n,h)), 



(54) 



where the "quasi-parity" p n ^ is any real or complex overlap between the eigenvectors 
before their biorthonormalization, 



A priori, the spectra of the eigenvalues of the operators H and/or A which satisfy 
the relations (49) may be either purely real or containing certain conjugate pairs of 
eigenvalues or containing certain "exceptional" points where the diagonalizability 
has been lost [11]. 

Whenever one contemplates the former option only, the reality of the eigenvalues 
opens a way towards a physical interpretation of the model where the eigenvalues of 
operators H and A may be observable and where the distribution of these eigenvalues 
can be experimentally measured. In the literature, a successful attempt of such a 
type has recently been completed for the spinless Klein-Gordon system [5]. In our 
present paper we made the first steps in the similar direction for spin one. For 
several alternative versions of the Hamiltonian H and associated helicity operator 
A we constructed certain "sufficiently simple" operators V playing the role of an 
indefinite metric which characterizes the relativistic Proca field of the massive vector 
bosons with spin one. 



((n, h\n, b!) = b~ n ,n'0~h,h'Pn,h- 
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5.2 Outlook 



Our present study has been inspired by the well known fact that in the massless 
limit of the vector bosons the construction of the physical metric becomes trivial. In 
particular, the operator of eq. (2) degenerates to a projector in a way explained by 
Gupta and Bleuler [26] (cf. also Appendix A below). In this sense our present study 
could be perceived as a preparatory clarification of several "PT— symmetry-related 
questions emerging during the canonical quantization of the massive model. 

Even in the absence of interactions we may expect that the construction of the 
positively definite 9 will be much less easy than in the electromagnetic massless 
case. Its details seem to represent an open problem up to now [18]. At the same 
time, the construction of 6 should be considered urgent at any spin. Its phenomeno- 
logical relevance stems from the necessity of a search for a compromise between the 
complexity of the field-theoretical calculations and the comparative simplicity of 
the various pragmatic semi- or non-relativistic models of the bound states (say, in 
a pionic atom, etc). We believe that the challenge of the construction of for spin 
one [25] might also re-open several new purely theoretical questions paralleling the 
ones encountered already in the zero-spin context of the Feshbach-Villars version of 
the first-quantized Klein-Gordon equation [5]. 
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Appendix A. Parallels between massive spin-one 
bosons and massless photons 

In a way resembling the description of the massless (i.e., electromagnetic) fields let 
us now change the notation and introduce an antisymmetric tensor F^ v such that 

F ij = — £ B 

F° j = -Ej, i,j = 1,2,3. 

In a fixed inertial frame equations (3) then acquire the form 



B = rot A, §T = - E -g r adA 

A = -m" 2 divE, |f = m 2 A + rotB. 



(55) 



An elimination of B and Aq gives 

dA ^ V 2 E q ^ dE 2A V 2 A g A 

where we abbreviated 

q = 2 grad div —V 2 . 

Finally, another abbreviation 

$ T = (u, v), u = E + imA, v = E-imA (57) 
converts equations (56) into the six-component evolution rule 

= & U = H TS 3>. (58) 

The matrix operator Hts is not Hermitian and fulfills the relation 

Hts = cr 3 H TS a3, ct 3 = ^ ^ 3 j 

shared by the spinless and spin-one cases [3, 27]. 



(59) 
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